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This  report  differs  greatly  from  the  work  which  I 
had  envisioned  at  the  onset.  The  plan  of  accomplishing 
a complete  performance  analysis  of  a K-ary  CPFSK  re- 
ceiver for  numerous  values  of  M and  n in  atmospheric 
noise  had  to  be  shelved.  The  main  reason  for  this  was 
the  large  amount  of  time  required  to  develop  an  ampli- 
tude probability  distribution  expression  to  approximate 
the  Crichlow  atmospheric  noise  model.  The  density 
function,  associated  with  this  expression,  was  required 
to  .evaluate  the  receiver  performance.  Therefore,  the 
receiver  was  evaluated  for  only  M=4  and  n=3. 

I would  like  to  express  my  sincere  gratitude  and 
thanks  to  my  patient  and  guiding  advisor.  Major  Carl 
of  the  Electrical  Engineering  faculty.  Despite  an 
extremely  busy  schedule  of  his  own,  he  always  found 
time  for  helpful  discussions.  The  discussions  and 
assistance  provided  by  Capt  Nejezchleb  were  also 
valuable  in  the  development  of  the  atmospheric  noise 
model  . 
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Cont inuous -phase  frequency  shift  keying  (CPFSK) 
is  potentially  an  attractive  modulation  scheme  with 
performance  better  than  Phase-shift  keying  (PSK)  in 
Gaussian  noise.  The  optimum  coherent  M-ary  CPFSK  re- 
ceiver has  been  developed  for  operation  in  Gaussian 
noise.  In  this  report  the  optimum  receiver  perfor- 
mance is  evaluated  assuming  operation  at  very  lov/ 
frequencies . 

Atmospheric  noise  is  the  most  prevalent  inter- 
ference when  operating  at  very  low  frequencies.  An 
analytic  amplitude  probability  distribution  (APD) 
function  is  developed  that  closely  approximates  the 
Crichlow  atmospheric  noise  APD  model  which  is  based 
on  empirically  measure d data.  The  analytic  expression 
for  the  APD  was  developed  by  combining  a Rayleigh  and 
a lognormal  density  and  obtaining  the  corresponding 
distribution  fi.T.nction.  By  correctly  choosing  values 
for  tv/o  parameters  in  the  distribution,  it  approxi- 
mates the  atmospheric  noise  APD  for  various  Vd  ratios 
( rms  vol  tage / ave rage  volr-age). 

The  performance  of  the  optimum  coherent  M-ary 
CPFSK  receiver  was  evaluated  for  M=4  operating  in  at- 
mospheric noise  with  specific  Vd  ratios  (4  & 5).  The 
evaluation  consists  of  finding  a probability  of  error 
bound  for  a sub-optimum  receiver,  which  is  valid  for 
high  si gnal - to-noi se  ratios  (SNR).  Results  indicate 


VI 


that  the  receiver  performance  in  atmospheric  noise 
becomes  increasingly  v/orse  compared  to  performance  in 
Gaussian  noise  as  the  SNR  increases. 


PERFORMANCE  ANALYSIS  OF  CONTINUOUS-PHASE 


FREQUENCY  SHIFT  KEYING  MODULATION 
IN  ATMOSPHERIC  NOISE 

I . Introduction 

There  is  a continuous  search  for  signal  modula- 
tion schemes  that  have  improved  performance.  Perfor- 
mance is  improved  if  a signal  can  be  transmitted  and 
received  with  the  same  error  rate  and  bandwidth 
requirements  but  at  a lower  si gnal - to-noi se  ratio 
(SNR).  Of  course,  the  source  and  type  of  signal 
interference  must  be  considered  in  evaluating  perfor- 
mance . 

Atmospheric  radio  noise  is  the  major  source  of 
errors-  when  operating  at  very  low  frequency  (VLF). 

The  physically  massive,  but  electrically  small,  trans 
mitting  antennas  required  for  operation  at  VLF  are 
also  a major  concern.  However,  VLF  communications 
systems  do  provide  a reliable  means  of  very  long 
range  surface-to-surface  communications.  Some  of  the 
attractive  charac  te  r i s t i c.s  of  radio  transmission  in 
the  VLF  band  are  (1)  very  low  propagating  signal 
attenuation  rates,  (2)  minimal  signal  fading,  and  (3) 
very  high  signal  phase  stability  (Ref  1:1).  These 
characteristic s make  VLF  communications  systems  a 
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good  choice  for  providing  command  and  control  communi- 


r 


j cations  in  times  of  war.  Therefore,  any  improvements 

i in  performance  of  these  systems  are  of  major  impor- 

tance . 

: Objectives  of  Research 

I 

i One  objective  of  this  research  is  to  develop  an 

t 

[ amplitude  probability  distribution  (APD)  function 

r 

I which  characterizes  atmospheric  radio  noise.  Crichlow 

f proposed  an  empirically  based  atmospheric  noise  dis- 

tribution model  based  on  attempts  to  identify  known 
distributions  which  could  be  used  to  represent  the 
probability  distribution  of  the  noise  envelope 
voltage.  The  distribution  function  developed  in  this 
report  should  provide  a model  that  corresponds  with 

i 

Crichlow *3  empiric  ally  based  model  (Refs  1:49-56;  | 

i 

3:5-8).  ,1 

The  second  objective  of  this  analysis  is  to 
evaluate  the  performance  of  continuous-phase  frequency  1 

shift  keying  (CPFSK)  modulation  in  the  presence  of 

■j 

atmospheric  radio  noise.  Performance  is  evaluated  by  I 

determining  the  symbol  error  probability  as  a function 
of  signal  power  and  noise  power. 

Background 

Although  there  have  been  many  years  of  experience 
in  VLF  radio  communications,  few  performance  analyses 
of  VLF  communications  systems  operating  in  atmospheric 

Ji 
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radio  noise  have  been  accomplished.  This  deficiency 
seems  to  exist  because: 

1.  Emphasis  was  placed  on  higher  frequencies  for 
radio  communications  where  atmospheric  radio 
noise  had  little  effect, 

2.  There  is  some  practical  difficulty  in  mathe- 
matically representing  bhe  atmospheric  radio 
noise  process,  which  is  non-Gaussian  and  non- 
stationary over  any  extended  period. 

3.  There  are  complications  involved  in  using 
nonlinear  receiver  techniques  to  optimize 
signal  demodulation.  Since  the  detailed 
statistical  characteristics  of  atmospheric 
noise  are  not  known,  optimum  receivers  were 
not  designed  for  operation  in  an  atmospheric 
noise  environment.  To  overcome  the  effects 
of  atmospheric  noise,  VLF  system  designers 
have  used  very  high  transmitter  powers  to 
compensate  for  inefficient  receivers  whose 
performance  is  drastically  degraded  by  impul- 
sive noise  (Ref  1:2-5). 

Recent  studies  indicate  that  CPFSK  modulation  can 
improve  the  efficiency  of  frequency  shift  keying  (FSK) 
communications  systems.  Since  FSK  is  gene>'ally  used 
in  VLF  communications  systems,  it  is  important  that 
any  modulation  scheme  that  improves  efficiency  be 
analyzed.  CPFSK  is  a modulation  scheme  wherein  the 
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phase  is  continuous  betv/een  frequency  shifts.  By 
keeping  the  phase  continuous,  transient  effects  are 
reduced  at  the  symbol  transitions,  and  this  provides 
spectral  bandwidth  advantages.  Another  important 
advantage  of  maintaining  a continuous  phase  is  that 
redundancy  is  imposed  upon  the  waveform.  Utilization 
of  this  reducdancy  allows  a decision  to  be  made  upon 
the  observation  of  several  symbols  rather  than  the 
more  common  approac  . of  making  a decision  by  looking 
at  each  symbol  independently  (Ref  4:644). 

Basi c Research  Strategy 

An  attempt  has  been  made  to  find  an  analytic  ex- 
pression for  the  atmospheric  noise  amplitude  density 
that  provides  results  comparable  with  the  Crichlow 
graphical  model.  Crichlow  proposed  that  the  APD  of 
atmospheric  noise  could  be  represented  by  two  straight 
lines  and  a circular  arc  when  plotted  on  probability 
paper  upon  which  the  Rayleigh  complementary  distribu- 
tion function  plots  as  a straight  line  (Refs  5:779; 
2:50;  6:714).  Wilson  has  developed  a density  function 
consisting  of  three  different  expressions  that  corres- 
pond to  different  regions  of  the  curve  (Ref  3:27-39). 
Although  this  density  function  consisting  of  one  ex- 
pression for  the  entire  curve  is  more  convenient  to 
use  in  determining  the  performance  of  various  re- 
ceivers. Research  has  been  based  on  the  idea  that 
the  ends  of  the  atmospheric  noise  APD  curve  approach 


two  specific  distribution  functions.  Therefore  the 
entire  APD  curve  may  possibly  be  represented  by  a 
distribution  function  that  is  based  on  a mixture  of 
density  functions  corresponding  to  distribution 
functions  approached  at  the  ends  of  the  curve. 

Once  an  atmospheric  noise  amplitude  density 
function  has  been  determined,  the  performance  of  a 
CPFSK  receiver  can  be  determined.  The  CPFSK  receiver 
is  optimized  for  M-ary  CPFSK  when  the  interference  is 
additive  white  Gaussian  noise  (Refs  7:1023-1036; 
4:644-652) . 

This  report  is  presented  by  chapters  in  the 
following  sequence.  Chapter  II  provides  a description 
of  the  coherent  CPFSK  modulation  scheme  and  of  the 
optimum  coherent  M-ary  CPFSK  detector  for  an  additive 
Gaussian  channel.  In  Chapter  III  an  atmospheric  noise 
model. is  developed  which  provides  results  that  closely 
approximate  the  Crichlov/  atmospheric  noise  model.  The 
performance  of  the  receiver  described  in  Chapter  II 
is  evaluated  in  Chapter  IV  assuming  that  the  interfer- 
ence is  atmospheric  noise.  The  atmospheric  noise 
model  developed  in  Chapter  III  is  used  in  this  evalua- 
tion. Chapter  V provides  conclusions  of  the  study  and 
recommends  areas  for  further  study. 
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II.  Coherent  M-ary  Continuous-Phase  Frequenc; 


Shift  Keying  Modulation 


The  objective  of  this  chapter  is  to  present  a 
description  of  M-ary  CPFSK  modulation  and  provide  a 
brief  description  of  the  optimal  coherent  M-ary  CPFSK 
detector  presented  by  Schonhoff  (Ref  4). 


Modulation  Description 

The  M-ary  CPFSK  received  signal  can  be  modeled 
1 1*1 

during  the  i^  signaling  interval  as: 


r(t)  = 1^1  cos  w^t 


dlTTh  ( t-  ( 1 -1 ) T 


( i-1 ) TctciT 


+ h ^ d . + 0 + 

j = l ^ 


n(  t ), 


where  E - energy  per  signaling  interval  (symbol 
period) 

u>^  - carrier  radian  frequency 

d^  - digital  information  signal  (±  1 , ±3  , . . . ±(M-1)) 
in  the  i^^  interval 

M - number  of  different  d^'s  (always  even) 
h - deviation  ratio 
0 - starting  phase 

n(t)  - additive  noise 

T - signaling  interval  length  or  symbol  period 


! 
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The  signal  is  constructed  so  that  the  separation 
between  two  adjacent  frequencies  in  the  M-ary  set  Is 
h/T  hertz.  For  coherent  CPFSK  the  starting  phase,  0, 
can  be  assumed  to  be  known  and  set  to  zero  with  no 
loss  of  generality.  The  shorthand  notation  developed 
by  Osborne  and  Luntz  (Ref  7:1024)  for  the  received 
signal  is  used  in  this  report. 

r(t)=S(t,d,D^)+n(t)  (2) 

where  d is  the  first  symbol  and  is  the  (n-1)  - tuple 

D = (d„,do,...d  ).  Equation  (1)  differs  from 
ot  ^ o n 

Equation  (2)  in  that  Equation  (2)  represents  the  re- 
ceived signal  over  n symbol  periods  and  Equation  (1) 
represents  the  signal  over  one  symbol  period.  Improved 
performance  can  be  obtained  with  CPFSK  modulation  by 
observing  the  received  signal  for  n symbol  periods  and 
making-  a decision  on  one  of  the  symbols  (Ref  4).  In 
coherent  detection  of  CPFSK  signals,  the  decision  is 
made  on  the  first  symbol. 

Cohe rent  Detection  o f M-ary  CPFSK 

The  receiver  design  uses  the  fact  that  the  carrier 
phase  during  the  i^  bit  time  depends  upon  the  data  in 
the  first  bit  time.  Schonhoff  (Ref  4:645-648)  devel- 
oped the  optimum  M-ary  coherent  detector  for  CPFSK 
assuming  that  the  interference  was  additive  white 
Gaussian  noise.  A block  diagram  of  this  receiver  is 
given  in  Figure  1.  In  this  figure  the  value  of  m 
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equals  which  is  the  number  of  possibilities 

for  D . The  receiver  correlates  the  received  waveform 

a 

with  each  of  the  m possible  transmitted  signals  begin- 
ning with  one  of  the  d^  values.  It  then  forms  a 
sum : 


m 

Lj^  = 5]^  exp(  2Xj^j/No)  , 


k = l 
J = 1 


• M 
. m 


(3) 


where  x,  . is  the  correlation  of  the  received  waveform 
kj 

with  the  si gna 1 wave  form.  The  decision  of  which 

data  value  was  sent  is  determined  by  the  largest  sum, 

S* 

The  evaluation  of  this  receiver's  performance, 
when  the  additive  noise,  n(t),  is  atmospheric  noise, 
requires  that  at  least  the  distribution  and/or  density- 
of  atmospheric  noise  be  known. 
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III.  Atmospheric  Noise  Model 


The  objective  of  this  chapter  is  to  develop  a 
mathematical  atmospheric  noise  APD  model  that  agrees 
with  the  Crichlow  model.  The  model  will  be  used  to 
determine  the  pe  r f o rmance  of  a CPFSK  receiver,  opti- 
mized for  an  additive  Gaussian  noise  channel,  while 
operating  in  an  atmospheric  noise  environment. 

Background 

Crichlow's  atmospheric  noise  APD  model  describes 
the  sta\.istical  properties  of  the  amplitude  distri- 
bution of  atmospheric  radio  noise.  The  data  used  to 
develop  this  model  consisted  of  measurements  of 
exceedance  probability 


rms 


(4) 


made  during  the  International  Geophysical  Year. 

Crichlow  noted  that  the  APD  could  be  approximated  by- 
two  straight  lines  and  a circular  arc  when  plotted  on 
paper  where  a Rayleigh  distribution  function  plots  as 
a straight  line  with  slope  equal  to  minus  one-half 
(Refs  5:778-780;  2:49-50).  An  example  plot  is  given 
in  Figure  2.  The  lower  portion  of  the  curve  represents 
low  voltages  with  high  exceedance  probabilities;  these 
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approach  a Rayleight  distribution.  The  upper  portion 
of  the  curve  represents  very  high  voltages  with  low 
exceedance  probabilities;  these  approach  a lognormal 
or  Weibull  distribution. 

The  APD  curve  varies  depending  upon  the  closeness 

of  atmospheric  noise'  generators  (lightning  discharges, 

etc).  As  the  lightning  discharges  occur  closer  to 

the  measuring  device,  the  atmospheric  noise  is  more 

impulsive.  This  noise  impulsiveness  is  described  by 

an  impulsiveness  form  factor,  , defined  to  be  the 

ratio  of  the  root  mean  square  (rms)  noise  envelope 

voltage  to  the  average  (ave)  noise  envelope  voltage 

(Ref  3:37-38).  Figure  3 shows  typical  APD  curves 

plotted  on  Rayleigh  paper  for  various  ratios. 

These  curves  show  that  as  the  V,  ratio  increases,  the 

d 

APD  of  atmospheric  noise  changes  in  a predictable  v;ay. 

Wilson  developed  the  expressions  for  the  three 
regions  suggested  by  Crichlow  and  determined  the 
boundary  points  between  the  regions  for  various 
ratios  (Ref  3:27-37).  He  represented  the  low  voltage, 
high  probability,  straight  line  portion  of  the  curve 
by  a Rayleigh  expression.  The  high  exceedance,  low- 
probability,  straight  line  portion  of  the  curve  was 
represented  by  a Weibull  expression.  The  center 
portion  of  the  curve  was  represented  by  a circular 
arc  expression.  From  these  three  expressions,  a 
density  function  was  developed  to  use  in  evaluating 
receiver  performance. 
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Using  Wilson's  noise  density  function  to  evaluate 


performance  of  various  receivers  is  not  convenient. 

It  is  advantageous  to  find  a single  expression  that 
represents  an  entire  APD  curve.  The  APD  curve  for 
each  ratio  could  be  obtained  by  changing  the  param- 
eters of  this  function. 

Rayleigh  — Lognormal  Mix  Distribution  Model 

Beckmann  noted,  as  Crichlow  and  Wilson  had,  that 
the  low  exceedance,  high  probability,  straight  line 
portions  of  Crichlov/'s  APD  curves  were  Rayleigh  dis- 
tributed. However,  he  found  that  the  high  exceedance, 
low  probability  portions  of  the  APD  lurves  could  be 
approximated  by  a lognormal  distribution  ( Ref  8:723-724 ) . 

Therefore  in  this  report  it  was  assumed  that  instan- 
taneous density  functions,  associated  with  both  the 
Rayleigh  and  lognormal  distribution  functions,  could  ! 

be  properly  combined  into  one  density  function  for  the  I 

entire  curve.  In  the  following  derivation  of  this 
mixed  density  function,  it  is  assumed  that  the  phase 
of  the  instantaneous  atmospheric  noise  voltage  is 
uniformly  distributed  and  independent  of  the  amplitude. 

This  assumption  is  made  for  expediency  since  there  is 

I 

no  known  relationship  between  noise  amplitude  and  j 

phase  . I 


If  the  atmospheric  noise  APD  is  Rayleigh,  then 
the  atmospheric  noise  amplitude  density  is  (Ref  9:49-50) 


(5) 


» 


P2.(z)  = -X  [-z^/A],  2 2-0 

where  z is  the  noise  amplitude  voltage.  If  the  instan- 
taneous atmospheric  noise  phase  is  considered  to  be 
uniformly  distributed,  the  density  function  of  the 
phase  is 

P0(0)  = H TTf  0 ?0  ? 2Tr  and  zero  elsewhere  (6) 

By  using  the  assumption  that  the  amplitude,  z,  and 
phase,  0,  of  atmospheric  noise  to  be  independent,  the 
joint  density  is  obtained  by  simply  multiplying  the 
densities,  (Equations  (5)  and  (6))  together. 

Pz,0(z*0)  = exp[-z^/A],  z ? 0 

0^0*2  (7) 

This  joint  density  function  can  be  transformed  into  a 
joint  density  function  with  in-phase  and  quadrature 
components  of  the  noise.  The  transformed  density  is 
sometimes  referred  to  as  an  instantaneous  density. 

The  in-phase  and  quadrature  components  of  this  density 
function  will  be  added  to  the  in-phase  and  quadrature 
components  of  an  instantaneous  density  function  ob- 
tained from  the  lognormal  distribution.  The  trans- 
formation is  accomplished  by  letting  x=z  cos  0 and 
^=z  sin  0 (Ref  9:10-12).  Therefore  the  instantaneous 
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density  is 


' |J|  Pj  0(  { f tan"^  (y/x)  ) 


where 


ItI  -i//  ^ 2 \ /^ 

|J  I = 1 / ( X + y j ^ 


the  Jacobian  of  the  transformation,  thus 


’x, = ST  <»'P  [- 


2 

X + 


(10) 


If  the  atmospheric  noise  APD  is  lognormal,  then 
the  atmospheric  noise  amplitude  density  is  (Ref  10:50) 


p^(z)  = 


.(2it)^ 


r_  ( ln(  z ) -a  ) ^ 

I 2s“ 


(11) 


Again  it  is  assumed  that  the  instantaneous  atmospheric 
noise  phase  is  uniformly  distributed  and  independent 
of  the  amplitude.  The  joint  density  is  obtained  by 
multiplying  densities,  Equati on  ( 6 ) an d Equation  (11), 


together  and  obtaining 


P2,0(Z,0) 


2sz  it(  2 it) 


exp  - 


r ( 1 n ( z ) - a ) 


0*022  TT 


(12) 
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This  joint  density  must  also  be  transformed  into  an 
instantaneous  density  function  as  was  done  in  the 
Rayleigh  case.  This  results  in 


5 5 V 

s(x  +y  )27r(27r)'2 


L 2s 


^2)^)-a)2- 


( 13  ) 


Assume  that  and  are  the  random  variables  in 
Equation  (10)  and  and  are  the  random  variables 
in  Equation  (13).  An  instantaneous  mixed  density 
function  with  random  variables  u and  v can  be  obtained 
by  letting  u = x^+ip  X = Z?  • Since  densities 

given  by  Equations  (10)  and  (13)  are  considered  inde- 
pendent, the  instantaneous  mixed  density  is  obtained 
by  convolving  the  densities  (Ref  11:188-189). 


Substituting  Equation  (13)  for  the  first  density  and 
Equation  (10)  for  the  second  density 


P.,  = 


21T  ( 2TT  ) '^s  ( x^  + y^  ) 


1 _r  (u-x)^  + 




( 1 n ( x^  + y ^ ) ^-a ) ^ 1 
L 2s^  ‘ 


( v-y ) 


dxdy 


(15) 


This  resulting  density  function,  Equation  (15), 
must  now  be  transformed  into  a Joint  density  of  magni- 
tude and  phase  random  variables.  This  is  done  by 
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letting  cos  0,  v=^  sin  0 and  applying  the  general 

transformation  equation  (Ref  12:96). 


p^  q(z,0)  =jjjp^  ^(z  cos  0,  z sin  0)  (16) 


where  jj|=  z.  the  Jacobian  of  the  transformation, 
Thus 


P_  ft(z,0)  = 

2Tr2(2TT)^sA 


2 2 
X + y 


;xp  - 


( ln(  + y^  ) '^-a  ) ^ 


•exp  - 


((z  cos  0-x)'^-t-(zsin9- 
A 


dxdy,  z » 0 


0 « 0 2ir 


(17) 


The  last  term  of  Equation  (17)  can  be  easily  expanded 


and  simplified  to 


r z^i  . v^i  fpz 

[-  — A 


(x  cos  0 + y sin  0)  (13) 


By  substituting  Equation  (18)  for  the  last  term  of 
Equation  (17)  the  Joint  density  becomes 


P„  c,(z.6)  = 5 — - — p — exp 

21T  “^  ( 2 IT) '^sA 


2 2 
X + y 


r (ln(x^  + y^)^-a)^  1 [ x^+  y^l  [2z,  . . 

' ^ ‘^2 ■ ^ expj^  — (x  cos  0-fysin0) 


dxdy,  z 


(19) 
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To  obtain  the  density  of  the  amplitude  of  atmos- 


pheric noise,  the  random  variable  0,  must  be  integrated 
out  of  density  Equation  (19). 


00 


The  last  term  in  Equation  (19)  can  be  simplified  to 
the  following  expression  by  a change  of  variable  and 
using  the  definition  of  a modified  Bessel  function. 


(22) 


Since  0 does  not  appear  in  Equation  (22),  it  is  easily- 
integrated  out  to  give 


3 


P (z)  = 

^ t \ /2 


2z 

-T7 exp 


(2  )'^sA 


[-  x] 


1 n ( r ) -a ) 


1- 


dr 


(23) 


The  root  mean  square  value  of  the  noise  amplitude, 
z,  that  is  lognormal  distributed,  is  (Ref  3:115-116). 


a s 

z = z z 
rms 


(24) 


Therefore,  if  the  root  mean  square  value  of  the  noise 

2 

amplitude  is  normalized  to  equal  1,  a = -s  . By  making 
this  substitution,  the  number  of  parameters  in  Equation 
(23)  can  be  reduced  from  three  to  two.  The  corres- 
ponding probability  of  exceedance  is 


Analytical  attempts  at  finding  a closed  form  solution 

to  Equation  (23)  and  (25)  failed.  If  the  substitutions 

2 2 

z=E,  z =E  „ , r=E  , s =a  , and  A=N  are  made  in 
rms  rms  o c 

Equation  (25),  the  resulting  equation  is  identical  to 
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Beckmann's  Equation  (5.19)  in  Reference  8.  This 
equation  is  the  final  solution  Beckmann  derived  to 
represent  the  APD  of  atmospheric  noise. 

To  determine  the  particular  values  of  parameters 
A and  s for  specific  ratios,  the  average  and  root 
mean  square  values  of  the  lioise  amplitude  are  solved 
using 

zp^ ( z ) dz  ( 26 ) 


z = square  root  of 
rms  ^ 


z p^ ( z ) dz 


(27) 


and  this  allowed  calculation  of  the  V,  ratio 

d 


V = 20  log 


rms 


ave 


(28) 


Equations  (25),  (26),  (27),  and  (28)  were  solved 

numerically,  using  Air  Force  Institute  of  Technology 
single  and  double  integration  programs.  These  two 
computer  programs  use  Simpson's  rule  and  are  identified 
as  SIMPS  and  SIMPD  respectively.  The  modified  Bessel 
function  of  the  first  kind  of  order  zero  (1^)  was 
evaluated  using  the  function  program  in  Appendix  A. 

This  program  is  developed  from  numerical  equations  for 
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evaluating  (x)  (Ref  13:377-378).  Numerical  solutions 
are  accurate  to  at  least  six  decimal  places. 

With  variables  r and  z integrated  from  0.001  to 
26.0,  the  probability  of  exceedance  is  plotted  for 
various  combinations  of  A and  s.  The  correct  combin- 
ation of  values  for  A and  s are  obtained  when  the 
following  two  requirements  are  met. 

1 ) 1 
rms 

2)  For  specific  ratios  the  APD  curves  obtained, 
using  Equation  (25),  must  match  the  APD  curves 
from  Crichlow's  model. 

Although  this  may  seem  relatively  simple,  it  is  highly 
improbable  that  the  exact  combination  of  values  will 
be  found  by  this  numerical  method.  Since  Equations 
(25),  (26),  and  (27)  were  not  numerically  integrated 

out  to  infinity,  the  most  that  can  be  hoped  for  is 
to  obtain  values  of  A and  s which  will  allow  close 
approximation  to  the  Crichlow  APD  curves.  Calculations 
for  the  higher  ratios  requires  a large  amount  of 
computer  processing  time,  thus  approximate  values  for 
A and  s are  obtained  only  for  ratios  of  four  and 
five.  These  values  are  listed  in  Table  I.  Since  the 
APD  curves  obtained  by  Wilson  appear  to  be  Identical 
to  the  Crichlow  curves,  Wilson's  curves  are  used  for 
comparison  purposes.  A computer  plot  comparison  of 
the  APD  curves  obtained  from  Wilson's  data  (Ref  3:45) 
and  from  Equation  (25)  are  given  in  Figures  4 and  5. 
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Figure  4.  Comparison  of  Distributions  Obtained  From 
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The  resulting  equation  for  the  probability  of  ex- 
ceedance is  thus  given  by 

00  00 

dT  z all  , 2ib/2)-l  f , 2,b/2l 

P zr(r  ) exp  - a(r  ) I- 

R o 


Again  no  closed  form  solutions  could  be  found  for 
Equations  (30)  or  (31).  Attempts  were  made  to  obtain 
the  proper  combination  of  values  for  a,  A,  and  b to 
obtain  approximations  for  the  various  APD  curves. 
However,  due  to  the  time  constraints  of  completing 
this  report,  the  many  possible  combinations  of  the 
parameter  values,  and  the  large  amount  of  processing 
time  required,  the  search  for  the  proper  combinations 
of  parameter  values  was  soon  dropped. 


J 

i 
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IV  . Pe rf o rmanc  e Evaluati on  o f the  Rece i ve  r 

The  objective  of  this  chapter  is  to  evaluate  the 
performance  of  the  receiver  described  in  Chapter  II 
when  the  signal  interference  is  caused  by  atmospheric 
noise.  The  1 ogno rmal -Ray le i gh  atmospheric  noise  mix 
model  developed  in  Chapter  III  is  used  in  this  evalua- 
tion. It  is  not  possible  to  analyze  the  performance 
of  the  receiver  exactly.  However,  bounds  on  perfor- 
mance can  be  determined  for  high  and  low  SNR.  Due  to 
time  constraints,  only  the  high  SNR  bound  is  developed 
in  this  report.  Developing  a low  SNR  bound  is  more 
complex  and  requires  the  development  of  a density 
K'  function  for  the  summer  outputs,  L's,  as  depicted  in 

Figure  1 . 

Hi gh-SNR  Cohe  rent  Bound 

Reference  Figure  1 the  M likelihood  parameters  can 
be  written  as  (Ref  4:646) 


For  large  SNR, 


where  x is  the  largest  of  the  x, 

— O "!■ 


Since  the  expo- 


nential function  is  a monotonic  function,  x is  an 

— o 

equivalent  parameter  to  investigate.  If  the  receiver 

simply  picks  the  largest  Xj.  • and  makes  r decision  on 

K J 

d,  it  is  a suboptimum  receiver. 


Density  o f Random  Vari ab 1 e Out  o f Cornel ators 

To  evaluate  the  performance  of  the  receiver,  the 
density  of  x v -i  be  determined. 

K j 


X,  . = 

-kj 


r( t)S( t ,k,D  )dt 
u 


(34) 


Assuming  that  the  transmitted  signal  is  S(t,d,D), 
where  D = {d2  , , . . . d^}  , Xj^j  is  given  by 


S(  t ,d,D)S(  t ,k,Dj  )dt  + / S(  t ,k,Dj  )n(t)  dt  (35) 


If  the  correlation  is  defined  to  be 


/IjdD  = ^ I S(t,d,D)S(t,k,D.)dt 


(36) 


Then  Equation  (35)  becomes 
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nT 


(37) 


= J"  S ( t , k , D j ) n(  t ) dt 


The  integral  in  Equation  (37)  is  now  defined  as  a 

random  variable,  £,  . 

K J 


-kj  ■ / S ( t , k , D J ) n( t ) dt 


( 38  ) 


The  density  of  c.  . is  equal  to  the  in-phase  or  quadra- 

K J 

ture  component  of  the  instantaneous  atmospheric  noise 
density,  since  Equation  (38)  can  be  considered  a 
narrowband  process.  The  instantaneous  atmospheric 
noise  density  is  given  by  Equation  (15).  By  letting 
r =x  + y , u=x,  and  v=y  Equation  (15)  is  simplified 


— ’ ^ A ITS  ( 2it  ) ^ 


00 

fl  r (In 

/-expj^ 


^ ^ 2.2 
( r ) + s ) 


(x^+y^-)-r^)|  J ^2r(x^  + y^)'^j 


dr  (39) 


The  density  of  £.  . is  then  just  the  marginal  density 

^ J 

derived  from  Equation  (39). 
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4 


Pc 

-kj 


*00 


( In ( r ) + 


2 . 2 
s ) 


2s‘ 


,2  2 2. 

( X +y  +r  ) 


^2r  ( X 


2 2. 

+y  ) 


drdy 


(40) 


This  density  has  zero  mean  and  is  symmetric  about  the 
origin.  The  density  of  x,  . is  given  by  (Ref  12:181) 

K J 


(41) 


In  words,  x,  . has  the  same  form  as  the  density  of  c 
-kj  - 

except  it  is  shifted  by 


Density  Evaluation . Equation  (41)  shows  that  the 

densities  p and  p are  identical  in  form,  but  they 
-dD  -kj 

have  different  means  which  are  determined  by  the 
correlation,  the  numerical  evaluation  of 

the  density  function.  Equation  (40),  r was  integrated 
out  to  15  and  y was  integrated  from  -10  to  +10  using 
the  Air  Force  Institute  of  Technology  double  integra- 
tion program,  SIMPD. 

The  density  function.  Equation  (40),  was  evaluated 
at  197  equally  spaced  points  between  -10  and  +10.  The 
limits  chosen  for  r,  y,  and  x are  picked  to  keep  inte- 
gration time  to  a minimum  and  still  have  at  least  two 
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M 


significant  figures  in  the  results.  Two  significant 
figures  are  considered  sufficient,  since  greater 
accuracy  cannot  be  obtained  from  Crichlow's  curves. 


Probability  of  Error  Derivation 


Assuming  that  d=v  and  D= { d2 , d^ » • • • } for  the 
transmitted  data  sequence,  the  probability  of  error 
is  given  by  the  following  equation. 


P[C‘|s(t,v,D)]  = P[any  Xj^j=-XvD^ 


(42  ) 


Applying  the  union  bound  (Ref  14:263-264). 


P[€'|s(t,v,D)] 


M m 

^ ^P[x,  .=»x  „] 

k?^v 


(43  ) 


M m 11 

P[C|  S(t,v,D)  I!  i / 

k=l  j=l  JJ 

k^v  -<*>  X 


-kj ’-vD 


( X , y ) dxdy 


(44) 


However,  the  joint  densities  of  different  combinations 
of  the  correlator  outputs  are  not  known.  Therefore,  a 
, worst  case  situation  will  be  assumed,  i.e.  the  Xj^^  are 
assumed  to  be  independent.  This  worst  case  condition 
will  give  an  upper  bound  on  the  probability  of  error 
since  the  correlation  between  the  Xj^ j ' s are  being 


ignored.  The  probability  of  error  is  then  given 

by 


P[f |S(t,v,D)  ]c 


M m 

Y.  Y 

k=l  j=l 

kj^v 


«00 


(y)dydx 

-kj 


(45) 


Then  by  averaging  the  conditional  probability  of  error 
over  all  possible  input  signals,  an  expression  for  the 
probability  of  error  is  obtained  that  is  independent 
of  the  input  signal. 


M m M 

Y Y Y 

k=l  j=l  v=l 
kj^v 


P„  ( y ) dy  dx  ( 46 ) 
-kj 


Numerical  Solution  of  Probability  of  Error 

The  bound  on  probability  of  error  for  M-ary  CPFSK 
is  solved  numerically  for  M = 4,  n = 3,  and  h = 0 . 8 . The 
value  used  for  h was  found  by  Schonhoff  to  be  the 
optimum  value  for  quaternary  modulation  and  n=3  in  an 
additive  Gaussian  noise  channel. 

Before  the  bound  on  the  probability  of  error  can 
be  determined,  Equation  (36)  for  the  correlation  must 
be  solved.  A symbo 1 -by-symbo 1 correlation  was  per- 
formed using  the  equations  developed  by  Schonhoff 
(Ref  4:647).  However  it  was  noted  that  Schonhoff 's 
equations  for  b and  d are  incorrect.  The  ir  h factor 
should  be  divided  out  of  each  equation,  making  them 
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n-1 

b = N + Y.  ^ ij 
1=2 


(47) 


n-1 

d = V + r d j 
£ =2 


(48) 


With  the  correlations  determined,  the  means  of 

densities  p (x)  and  p (x)  are 
-vD  -kj 


”SvD  ’ 

j = D 


(49) 


k 


■ /ijdD"':’  >'  * 1 


(50) 


f.  i 


To  obtain  a numerical  solution  to  the  inequality 
Equation  (46),  the  integrals  are  changed  to  summations 
with  limits  which  will  cover  the  non-zero  area. 


M m M m 

ptn-  T H Z L 

k=l  j=l  v=l  D=1 

k9^v 

M +10  +10 

^kj 

x=M„  -10  — vD  y=x  — kj 


-vD 


(51) 
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Since  the  numerical  solution  of  inequality 
Equation  (51)  requires  a large  amount  of  computer 
time,  the  probability  of  error  bound  is  calculated 
by  averaging  over  8 randomly  picked  input  signals  in- 
stead of  the  64  total  possible  signals.  Probability 
of  error  solutions  obtained  by  averaging  over  8 
randomly  picked  input  signals  were  compared  with 
solutions  obtained  by  averaging  over  different  groups 
of  8 randomly  picked  input  signals.  This  comparison 
revealed  that  the  probability  of  error  differed  only 
in  the  third  significant  digit.  Therefore  it  is 
assumed  that  this  approximation  does  not  effect  the 
results  appreciably.  The  inequality  Equation  (51)  is 
reduced  to 


M m 


prri«  V V ^/Over  8 arbitrary^  ^ 
i i..Mnput  signals  ^ 


M +10 


P„  (x) 

^-10 
— vD 


M +10 

L Pv  (y)AyAx 

y=x  ^kj 


(52) 


Plotti ng  Receiver  Performance 

Receiver  performance  is  displayed  by  plotting 
probability  of  error  bound  versus  SNR.  The  SNR  is 
expressed  as  the  energy  per  symbol  interval  divided 


by  the  variance  of  the  atmospheric  noise  density, 
Equation  (40).  It  is  assumed  that  the  energy  in  each 
symbol  interval  was  normalized  to  one  and  therefore 
there  is  equal  energy  in  each  symbol  interval. 

Performance  plots  for  V,  = 4.0  and  V,  = 5.0  are 

d d 

shown  in  Figures  6 and  7,  respectively.  The  perfor- 
mance is  also  plotted  for  Gaussian  noise  with  variance 
equal  to  the  variance  of  the  respective  noise  densities 
used  for  comparison. 

Figures  6 and  7 clearly  indicate  that  the  upper 
bounds  on  probability  of  error  for  atmospheric  noise 
and  Gaussian  noise  come  closer  together  as  the  SNR 
decreases.  At  SNR  = 0 the  upper  bound  for  each  is 
almost  identical.  The  implication  of  the  performance 
curves  is  that  if  the  receiver  is  used  in  an  additive 
atmospheric  noise  environment,  then  a higher  SNR  will 
have  to  be  mciintained  to  obtain  the  same  probability 
of  error  as  predicted  for  its  use  in  a Gaussian 
channel . 

It  is  noted  that  the  Gaussian  noise  curves  for 
M = 4,  n=:3,  and  h = 0.8  are  identical  for  different 
variances.  This  result  is  expected  since  the  proba- 
bility of  error  is  plotted  versus  SNR.  However,  the 
Gaussian  noise  curves  obtained  in  Figures  6 and  7 are 
not  identical  to  Schonhoff's  curve  (Ref  4:648).  A 
curve  comparison  is  given  ir;  Figure  8,  which  shows 
that  for  a SNR<2.0  Schonhoff  obtains  a lower  probability 
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Figure  6.  Probability  of  Symbol  Error  for  Quaternary 


Coherent  CPFSK, 


= 4.0 
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Figure  7.  Probability  of  Symbol  Error  for  Quaternary 
Coherent  CPFSK,  V.  = 5.0. 
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Figure  8.  Probability  of  Symbol  Error  Comparison 
With  Schonhoff's  Results. 
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V . Conclusions  and  Recommendations 

Conclusions  and  recommendations  cover  two  separate 
areas;  the  atmospheric  noise  model  and  the  coherent 
M-ary  CPFSK  rece 1 ve r pe rf ormance . 

Atmospheric  Noise  Model 

The  Rayleigh-lognormal,  mixed  distribution  model 
closely  represents  the  Crlchlow  APD  model  of  atmos- 
pheric noise  If  the  correct  parame-te^  values  for  A 
and  s are  found  for  each  particular  ratio.  The 
present  difficulty  in  using  this  model  is  the  trial 
and  error  method  used  to  obtain  the  correct  parameters 
for  various  ratios. 

Further  research  should  be  done  to  identify  a 
relationship  between  the  ratio  and  the  parameters 
A and  s so  a parameter  set  could  be  easily  found  for 
any  ratio.  Research  also  needs  to  be  continued  to 
find  an  atmospheric  noise  distribution  model  which  has 
a closed  form  solution. 

Cohe rent  M-ary  CPFSK  Receiver  Performance 

As  seen  in  Figure  6 and  Figure  7,  the  probability 
of  error  for  this  receiver  is  always  greater  when 
operating  in  an  atmospheric  noise  environment.  Wilson's 
evaluation  of  phase  shift  keying  indicated  that  at  low 
SNR,  a PSK  receiver  may  perform  better  in  atmospheric 


noise  than  in  Gaussian  noise  (Ref  3:100-105). 


Although  results  in  this  report  do  not  fully  justify 
the  same  conclusion,  it  still  may  be  true.  Performance 
results  of  a M-ary  CPFSK  receiver  at  low  SNR  may 
resemble  Wilson's  results  if  a performance  bound  for 
SNR  is  developed. 

There  are  many  areas  that  require  further  re- 
search in  evaluating  M-ary  CPFSK  receiver  performance. 
Five  major  areas  are: 

1.  Need  to  obtain  a performance  bound  for  low  SNR. 

2.  Neet  to  evaluate  receiver  performance  for 
different  values  of  M and  n. 

3.  Need  to  find  the  optimum  value  of  h (deviation 
ratio)  for  different  values  of  M and  n in 
M-ary  CPFSK  when  the  interference  is  atmos- 
pheric noise. 

4.  Need  to  analyze  the  performance  of  a non- 
coherent M-ary  CPFSK  receiver. 

5.  Need  to  develop  optimum  coherent  and  noncoher- 
ent M-ary  CPFSK  receiver  for  atmospheric  noise. 
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The  computer  program  used  to  evaluate  the  modified 
Bessel  function  of  the  first  kind  order  zero  is  given 
below. 


C FIO=MCDIFIED  BESSEL  FUNCTION  OF  FIRST  KINO  ORDER  0 
FUNCTION  FTO(X) 

DIMENSION  XL  (6  ) ,XH(8) 

DOUBLE  W 

OITA  XL /3. 5 158  E29. 3. 08R'54?  4,  1.2Q57i|9?,  0.26  5973?,  0.0  360  768, 
10. 0 045813/, XH/ 0.0 13? 3592, 0.00 225315, -O.CD 157555, 0.03916281, 
2-0.02057706,0. 02635537,-0.01647633, 0.00392377/ 

W = X 

yA  = AR^  f y ) 

IF  (XA.GT. 0.001)  GO  TO  10 

FIO=1.0 

RETURN 

10  IF  (X.GT.O)  GO  TO  20 
X=XA 

21  T=X/3.75 

IF  (T.6T.1.0)  GO  TO  40 

Fro=i.o 

DO  30  U=l,6 

FIO=FI  0+XL  (J)*  (T**(2'^J)) 

30  continue 

RETURN 

40  FI0=0. 39894228 

TI  = 1 .0/T 
DO  50  J=l,8 
FIO=FIO+XH(J)*  (TI-»*J) 

50  CONTINUE 

FIO=FIO*OEXP(W)» (l.O/SQRT (X) ) 

RETURN 

END 
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